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ABSTRACT 



Context. Classical Be stars are hot non-supergiant stars surrounded by a gaseous circumstellar disk that is responsible for the observed 
ON ' infrared (IR) excess and emission lines. The influence of binarity on these phenomena remains controversial. 

Aims. We followed the evolution of the environment surrounding the binary Be star S Scorpii one year before and one year after the 
201 1 periastron to check for any evidence of a strong interaction between its companion and the primary circumstellar disk. 
Methods. We used the VLTI/AMBER spectro-interferometric instrument operating in the K band in high (12000) spectral resolution 
to obtain information on both the disk geometry and kinematics. Observations were carried out in two emission lines: Bry (2. 172 yum) 
and He I (2.056 ^m). 

Results. We detected some important changes in S Scorpii's circumstellar disk geometry between the first observation made in April 
2010 and the new observation made in June 2012. During the last two years the disk has grown at a mean velocity of 0.2 km s . This 
is compatible with the expansion velocity previously found during the 2001-2007 period. The disk was also found to be asymmetric 
5— I , at both epochs, but with a different morphology in 2010 and 2012. 

Conclusions. Considering the available spectroscopic data showing that the main changes in the emission-line profiles occurred 
quickly during the periastron, it is probable that the differences between the 2010 and 2012 disk geometry seen in our interferometric 
data stem from a disk perturbation caused by the companion tidal effects. However, taking into account that no significant changes 
have occurred in the disk since the end of the 201 1 observing season, it is difficult to understand how this induced inhomogeneity has 
been "frozen" in the disk for such a long period. 

> ■ 
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00 ' 

1 . Introduction boundary (photosphere/disk interface) rotating close to the criti- 

■ cal velocity (V c ). The expansion was found to be approximately 
, We initiated an interferometric follow-up of the binary Be 2kms _1 

— , ' star 6 Scorpii which has been regularly observed with 

• • ; the VLTI/AMBER (Petrov et al. 2007) and CHARA/VEGA ,„ . . , . . f „ c , m , „ n . 

L> /l/r j . , »nnm ■ » » inm t i\/r n a » With the measured vsim of 175 kms (Berbacca & Pennotto 

. — ■ (Mourard et al. 2009) instruments since 2007. In Meilland et , , ...... , „ „_ „. 

i /omi\ u n t a * a s c ••> 1970), and the determined inclination angle of 28 ±8 we m- 

al. (2011), hereafter Paper I, we managed to refine 5 Scorpii s _ , & , „. . . , 

binary orbit and the orbital period was the only orbital element fe "" ed tha " he star rotates between 60% and l05 J° of f cntlcal 
3 , • / 'ii .i A -Ar t e tU i . • j r velocity. Moreover, taking into account a possible underestima- 
te ■ which was significantly different from the one determined by / , ... 6 . . . . /rr , 

rp . , /o AAri N ^ t1 j . ■ . tion of the vsim because of gravity darkening (Townsend et al. 

Tango et al. (2009). Consequently, we determined that the next , , ... , „ . , , . 

periastron should take place around 201 1 July 5 (±4 days). This ™) the lower limit on the stellar rotation would rise to about 

, C J\ j u m . , /omn , . 70 % V c . Finally, we measured the extension in the Ha line, i.e., 

result was also confirmed by Tycner et al. (2011). Thanks to . _ , e , , , , . . 

.• eu- u \- 1 a i i 4.8 ± 1.5 mas, which is about the same as the binary separation at 

the combination of high spectral and spatial resolution we con- , . ' . _ , , 

strained the geometty and the kinematics of 8 Scorpii's circum- the Penasfron, i.e., 5.9 mas. Thus, we suggested that a strong in- 

. „ • . c i .1 . .1 i • • • • • teraction between the already ejected matter and the companion 

stellar environment. We confirmed that the line emission ongi- ... , . , . 

. e A . • i j* i a n a ■ t au i ti should occur during the periastron passage, 

nates from a dense equatorial disk fully dominated by rotation. & f f & 

We found that the rotation appears to be Keplerian, with an inner 

This letter is organized as follows: in Sect. 2 we briefly summa- 

Send offprint requests to: ame@oca.eu rize other recent observations of 6 Scorpii; in Sect. 3 we present 

* Based on observations made with ESO telecopes at La Silla Paranal our VLTI/AMBER observations and the data reduction. In Sect. 

Observatoiy under GTO programme IDs 089.D-0105(A) and 089.D- 4, the data are analyzed and compared to those acquired in April 

0105(B) 2010, one year before the periastron, and finally, a short conclu- 

** Fig. [3]is only available in electronic form via www.edpsciences.org sion is drawn in Sect. 5. 
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2. Recent observations of 6 Scorpi 



Table 1. VLTI/AMBER log for the 2012 6 Scorpii observations. 



the star S Scorpii (7 Scorpii, HR 5953, HD 143275) was ob- 
served by Millan-Gabet et al. (2010) using the CHARA/MIRC 
instrument in the H band. They resolved an elongated disk with a 
Gaussian FWHM of 1.18 x 0.91 mas. In the K band the disk was 
only marginally spatially resolved and they estimated an uncer- 
tain continuum disk FWHM of 0.7 + 0.3 mas. When taking into 
account the continuum/line flux ratio they obtained much larger 
sizes for the line emission regions, i.e., 2.2 + 0.4 mas for He I and 
1.9 + 0.3 mas for Bry. In Paper I, we resolved the circumstellar 
disk in the Ha (FWHM=4.8 ± 1.5 mas), Bry (FWHM=2.9±0.5 
mas), and the 2.06 /mi He : (FWHM= 2.4 + 0.3 mas) lines, as 
well as in the K-band continuum (FWHM ~ 2.4 mas). 

In a recent paper Che et al. (2012) imaged the disk distortion 
near the periastron. They revised the binary orbit and studied 
the primary disk properties during the periastron with NPOI and 
CHARA/MIRC interferometers. They were able to fit a global 
model where the secondary follows a Keplerian orbit, as already 
mentioned in Paper I, and the disk properties were found to be 
stable through all their observations. The estimated the new peri- 
astron passage to be UT 201 1 July 3 07:00 + 4:30 which is within 
the error bar of the prediction made in Paper I. They also esti- 
mated the mass of the secondary, based on the revised binary or- 
bit and radial velocity (RV) measurements from Miroshnichenko 
et al. (200 1 ), to be ~ 6 Mq with the primary mass 13.9 Mq esti- 
mated from v sin i, apparent Teff, and V-band photometry mea- 
surements. The mutual angles between the disk mid-plane and 
orbital plane were estimated to be either ~119°or~171°. From 
the fitted global model, they also found that ~ 3% of the H-band 
flux comes from a fully resolved envelope. The primary disk 
was found to be mainly symmetric and stable, contributing to 
71.4 + 2.7% of the total H-band flux. This implies a quiescent 
inner disk and no ongoing material outflow after the periastron 
which is very different from what we obtained in the K band. 

On the other hand, Rivinius et al. (2012) conducted a spec- 
troscopic study of 5 Scorpii emission lines during the periastron. 
They found that Ha and H/3 profiles were strongly modified at 
the periastron. The major modification was that the Violet/Red 
(V/R) ratio went from V/R>1 to V/R<1 in only a few days. They 
argued that this was due to the direct tidal effect of the compan- 
ion on the primary disk. 

We note that we looked at many Ho- and H/3 observations 
made by amateurs available in the BeSS database^ and that the 
line profiles did not change significantly from the state detected 
by Rivinius et al. (2012) after the periastron. 
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Fig. 1. 6 Scorpii (u,v) 491an coverage obtained in 2012 with 
VLTI/AMBER observations centered on Bry (left) and 2.06 /im 
He I (right) emission lines. 



Date Time 




L. (m) 


P.A.(°) 


Observations centered at 2.17 yum (Bry) 


2012-06-22 23:46 
2012-06-23 00:32 
2012-06-23 01:21 
2012-06-23 02:05 


D0-H0-G1 
D0-H0-G1 
D0-I1-H0 
D0-I1-G1 


53/71/ 63 
57/71/ 67 
79/ 40/ 62 
82/46/71 


48/ 170/ 124 
57/ 175/ 126 
94/ -36/ 63 
98/- 143/ 132 


Observations centered at 2.06 //m (He/) 


2012-06-23 23:16 
2012-06-24 00:07 
2012-06-24 00:52 
2012-06-24 01:35 
2012-06-24 02:18 


D0-H0-G1 
D0-I1-H0 
H0-I1-G1 
D0-I1-G1 
D0-I1-G1 


50/71/ 60 
70/ 39/ 56 
40/45/71 
80/46/71 
82/47/71 


43/ 167/ 123 
87/ -40/ 53 
142/-152/ 177 
95/- 146/ 130 
99/-141/ 134 



3. Observations and data reduction 

The star 6 Scorpii was observed in June 2012 during two 
half-nights with the VLTI 1.6m auxiliary telescopes (AT) us- 
ing the AMBER beam recombiner (Petrov et al. 2007). The 
observations were carried out using high spectral resolution 
mode (R= 12000) centered on either Bry or Hei (2.06 yt/m) for 
each of the two half-nights. In order to enhance the signal- 
over-noise ratio (S/N) by enabling longer exposure time, the 
fringe tracker FINITO was also used. The star HD 139663 
(D = 2.133 + 0.147 mas, Lafrasse et al. 2010) was used as an in- 
terferometric calibrator. 

The data were reduced using the standard AMBER data 
reduction software amdlib v3.0.3bl (see Tatulli et al. 2007 
and Chelli et al. 2009 for detailed information on the AMBER 
data reduction). The average-raw complex visibility and closure 
phase was determined using the standard method, i.e., keeping 
the 20 % of the frames with the higher S/N ratio. The interfer- 
ometric calibration was then done using custom-made scripts^ 
described in Millour et al. (2008). The (u,v) plan coverage for 
all the observations is plotted in Fig.Q] 

4. The disk geometry and kinematics 

4.1. Qualitative comparison of the 2010 and 2012 data 

The Bry and Hei absolute visibility show a strong variability 
with an amplitude of about 0.2. Assuming that the binary sepa- 
ration at the epoch of the observations should be 90.8 mas with 
a position angle (PA) of 156.08°, according to Paper I orbital 
parameters, these variations are compatible both in terms of am- 
plitude and frequency with the binary modulation. However, the 
accuracy of the absolute visibility of VLTI/AMBER data taken 
in HR mode is not sufficient to significantly refine the orbital 
parameters. 

The differential visibilities and phases, closure phase, and 
line profile from the 2010 Bry observation are plotted in Fig. [2] 
In the same figure, one single Bry 2012 observation with similar 
projected baselines is overplotted. The full 2012 dataset for Bry 
and He I emission lines is presented in Fig.[3](online). 

The visibility and phase variations shown in Fig. [2] are typ- 
ical for a rotating disk. As discussed in detail in Meilland et al. 
(2012), the morphology of visibility variations for a purely ro- 
tating disk depends on the baseline orientation. Considering a 
non-fully resolving baseline, the variation is V-shaped for a base- 
line aligned with the disk's minor axis (baseline 1 of Fig |2]l and 
W-shaped for a baseline aligned with the major axis (baseline 



available at https : //www . oca . eu/spip . php?article 562 
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Fig. 2. 2010 (in red) and 2012 (in blue) Bry VLTI/ AMBER data for a similar projected triplet of baselines. 



2 of Fig |2). Baselines with intermediate orientation exhibit in- 
termediate morphologies (see, for example, the weak W-shaped 
variation for the third baseline shown in Fig. |2). 

The S-shaped phase variations also strongly depend on the 
baseline orientation. In the case of a purely rotating disk, the 
amplitude is maximum for a baseline aligned with the major 
axis and zero for a perpendicularbaseline. Thus, from the simple 
analysis of data shown in Fig.|2]it seems that the disk orientation 
did not significantly change between 2010 and 2012 with the 
major axis being oriented roughly in the direction North-South. 

However, Fig. [2] also clearly shows important changes in 
the visibilities, phases, and line profile between 2010 and 2012. 
Some of them favor the hypothesis of a significant growth of the 
equatorial disk between the two epochs: 

• The visibility drop is clearly deeper in 2012. 

• The 2012 S-shaped phase variations are larger. 

• The 2012 line profile equivalent width (EW) is larger. 

• the double-peak separation is narrower in 2012. 

Nevertheless, other major changes seen in Fig. |2]are not re- 
lated to a growth of the disk: 

• The V/R ratio went from 1.2 in 2010 to 0.9 in 2012. 

• The W-shaped-visibilities peak ratio is reversed. 

• The 2012 S-shaped phase variations are asymmetric. 

Such changes of asymmetries in the spectrum, visibilities, 
and phases across the Bry line, are clear evidence of changes 
in the circumstellar matter distribution around the star. In the 
case of a rotating disk, these asymmetries usually stems from the 
inhomogeneity of the disk. The origin and nature of the disk's 
inhomogeneities and their evolution around the periastron are 
discussed in Sect. 4.3. 

We note that all this evidence of the evolution of the disk 
around 6 Scorpii between 2010 and 2012 are seen not only in 
the data plotted in Fig. [2] but also on all the Bry and He I dataset 
(see Fig.[3]online). 

4.2. Modeling of the 2012 data using the kinematic model 

To model the wavelength dependence of the visibility, differen- 
tial phase, and closure phase in the observed emission lines we 
used the simple kinematic model developed for fast fitting of an 
expanding and/or rotating thin equatorial disk already described 
in paper I. The model is presented in more detail in Delaa et al. 
(201 1) and Meilland et al. (2012). 



We used the best- fit model parameters of the 2010 data from 
Paper I as a starting model. As expected, we obtained a very 
unsatisfactory fit with a x 2 r of 24.7. We then assumed that the 
stellar and kinematics parameters should not change between 
2010 and 2012, and that only the Bry and Hei disk exten- 
sion and line equivalent width (EW) should have evolved. In 
this case we found a best-fit model for FWHM Bry = 8.0 + 0.5 R* 
(4.2 + 0.3 mas), EW Br? = 7.2 + 0.3 A, FWHM/m = 5.2 + 0.3 R* 
(2.8 + 0.2 mas), and EW He! = 9.1 ±0.3 A. However, for this set 
of values the^ 1S still high, i.e., 17.4. 

We then performed a fit after releasing constraints on the 
stellar inclination angle (i), disk major axis orientation (PA), 
and disk kinematics, i.e., rotational velocity at the photosphere 
(V rof ) and exponent of the rotational law in the disk (/?). Thus, 
the only fixed parameters for this set of models are the distance 
taken from van Leeuwen (2007) and the stellar radius taken from 
Paper I. The best-fit model parameters are presented in Table [2] 
and the model visibility, phases, and line profiles are overplotted 
on Fig. [3] (online only). 

The^ for this best-fit model, however, is still high, i.e., 14.6. 
By separating the fit on each line we found that the He I data 
have a significantly better fit (i.e., ;^ = 8.6) than the Bry data 
Otfr =22.7). For both lines, the largest deviations between the 
observations and the data comes from the asymmetrical shape of 
the line profile as well as some visibility and phase variations. 
This effect will be discussed in the next section. 



Table 2. The 2010 and 2012 best-fit kinematic models. 



Parameter 


(unit) 2010 


2012 




Global geometric parameters 




R* (fixed) 


(Re) 8.5 


8.5 


d (fixed) 


(pc) 150 


150 


i 


(deg) 30 ± 8 


26 + 5 


PA 


(deg) -12 ±7 


-14 + 4 




Global kinematic parameters 




v rot 


(km.s- 1 ) 500 ± 50 


500 ± 50 


P 


0.5±0.1 


0.5 + 0.05 




Bry disk geometry 




FWHM Bry 


(R*) ' 5.5 ±1 


7.8 + 0.3 


EW Bry 


(A) 6.5+0.5 


6.2 + 0.2 




He; disk geometry 




FWHMhci 


(R*) 4.5+0.5 


4.9 + 0.2 


EWhci 


(A) 8.5+0.5 


8.8 + 0.3 


xi 


6.2 


14.6 
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We also note a second source of deviation. Our best-fit model 
clearly underestimates the visibility for baselines close to the 
disk's minor axis. However, using a larger disk extension, we 
would overestimate the visibility for the other baselines and the 
over-resolution of the disk would strongly modify the shape of 
the phase variations. Such underestimation of the minor-axis vis- 
ibility was already detected for a Col, one of the Be stars ob- 
served with VLTI/AMBER (Meilland et al. 2012). However, its 
origins remain unknown. 

4.3. The disk asymmetry 

Okazaki (1997) proposed that the cyclic variation of the V/R 
peak ratio with a period of several years originates in the pre- 
cession of one-armed density waves in the disk. This assump- 
tion was compatible with spectro-interferometric observations 
of y Cas (Berio et al. 1999), K CMa (Meilland et al. 2007), 
and £ Tau (Caricofi et al. 2009). This could also explain the 
change of the peaks ratio of the W-shaped visibility seen in 
Fig. |2] However, according to observations made by Rivinus et 
al. (2012), the V/R changes for 5 Scorpii occurred during a very 
short timescale of a few weeks around the periastron. Thus, it 
seems that the hypothesis of the one-arm oscillation cannot ex- 
plain the current disk properties. 

We also note that, using our simple axisymmetric model, 
we obtained a better fit on the He I data than on the Bry data. 
Considering that the Bry emission line is formed further from the 
primary, it is supposed to be more perturbed by the companion 
tidal forces than the He I line. Consequently, the worse fit of the 
Bry data might be a clue of a more important tidal perturbation. 
One can also argue that the broad-band measurements taken by 
Che et al. (2012) could not detect any tidal effect as the H-band 
continuum stems from an even smaller part of the disk than the 
He I line. Nevertheless, in the frame of an asymmetry induced by 
tidal effects in June-July 2011, it is also difficult to understand 
why the V/R remained almost constant between July 2011 and 
our observations in June 2012. One should expect the induced 
inhomogeneity to be either in co-rotation or precession, causing 
cyclical V/R variations with timescales ranging from weeks to 
years that are not seen in the available spectroscopic data. 

4.4. Evolution of the disk extension 

Our VLTI/AMBER data clearly gives evidence that the Bry 
emission FWHM increased by about 40% (i.e., 2.3 + 1.3 R*) be- 
tween 2010 and 2012. Assuming that this was directly cause by 
a physical growth of the circumstellar disk we derived a mean 
expansion velocity of 0.20+0.1 1 km s _1 . This value is similar to 
the one previously determined for 6 Scorpii by Miroshnichenko 
et al. (2003) for the 2001-2003 period, i.e., 0.4 km s -1 , and the 
one derived in Paper I for the 2000-2005 and 2005-2007 peri- 
ods, 0.24 km s~' and 0.19kms _1 , respectively. It is also com- 
patible with the one determined by Kanaan et al. (2008) for 
the transient disk surrounding the binary Be star Achernar, i.e., 
0.27 ±0.08 km s- 1 . 

The He I line, however, did not exhibit the same variation as 
the Bry line. Its growth, of about 8%, remains within the un- 
certainties of the derived extension. One possible explanation 
for this difference is that He I has a higher excitation energy. 
Consequently, it is produced by denser and hotter material than 
Bry. Thus, even if the disk is increasing in size as a function of 
time the He I emitting region should remain compact, whereas 



the Bry may be produced in the newly formed outer part of the 
disk. 

In a recent paper, Stee et al. (2012) measured the disk exten- 
sion in these two lines for the relatively steady disk of the Be 
star HD 1 10432. The values derived for their measurement using 
the same methods showed that the He I emission is about 25% 
smaller than the Bry emission. This value is between the one we 
found for 6 Scorpii for our 2010 (18%) and 2012 (37%) models. 

5. Conclusion 

Using the unique high spectral and spatial resolutions of 
the VLTI/AMBER instrument, we studied the evolution of 
5 Scorpii's circumstellar environment one year before and one 
year after the 201 1 periastron passage. At both epochs, the data 
clearly shows that the ejected matter is located in a thin equato- 
rial disk in Keplerian rotation as it is for other observed Be stars. 
Since 2000, the disk seems to have grown with a quasi-constant 
velocity of about 0.2kms _1 and was still growing during the 
2010-2012 period. Such speed is compatible with the viscous 
excretion disk model proposed by Lee et al. (1991). 

We found evidence that the disk was asymmetric at both 
epochs of our observations, but with a different morphology in 
2010 and 2012. According to spectroscopic observations made 
by Rivinius et al. (2012), the major modifications occurred near 
the periastron. Thus, the observed asymmetry cannot originate 
from a one-arm oscillation as proposed by Okazaki (1997), but 
must stem from tidal effects induced by 6 Scorpii's companion. 
However, according to available spectroscopic data and our 2012 
observations, the disk asymmetry seems to be "frozen" since the 
periastron. Finally, long-term spectroscopic and interferometric 
follow-up are needed to understand both the nature of the inho- 
mogeneity after the periastron and its link with the tidal effects 
of the companion. 
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Fig. 3. The fuU dataset for <5 Sorpii 2012 VLTI/ AMBER observations. The red lines represents the VLTI/AMBER data, and thi 
green ones, the best-fit symetric kinematic model. 



